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Abstract 

The purpose of the present paper is two-fold. On the one hand, we review some recent studies on 
the low-temperature strong-held thermodynamic properties of frustrated quantum spin antiferro- 
magnets which admit the so-called localized-magnon eigenstates. One the other hand, we provide 
some complementary new results. We focus on the linear independence of the localized-magnon 
states, the estimation of their degeneracy with the help of auxiliary classical lattice-gas models and 
the analysis of the contribution of these states to thermodynamics. 

PACS numbers: 75.10.Jm; 75.45. +j; 75.50.Ee 

Keywords: frustrated Heisenberg antiferromagnets, strong magnetic field, localized-magnon states, classical 
lattice-gas models 
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I. INTRODUCTION 



Antiferromagnetically interacting quantum Heisenberg spins on regular geometrically 



ar geometric 
US 3- On 



the 



frustrated lattices have attracted much interest over the past few decades 
one hand, there has been tremendous recent progress in synthesizing corresponding magnetic 
materials 4j]. On the other hand, geometrically frustrated quantum antiferromagnets are 
an excellent play-ground for studying novel quantum many-body phenomena. We mention 
quantum spin-liquid phases, order-by- disorder phenomena, lattice instabilities to name just 
a few. Application of an external magnetic field to a frustrated quantum Heisenberg antifer- 
romagnet introduces a new competition between interactions in the spin system that may 
lead to further interesting phenomena. As an example we mention the half-magnetization 
plateau stabilized by structural distortion in the pyrochlore lattice 

The theoretical investigation of frustration effects in quantum spin antiferromagnets 
usually meets new difficulties; e.g., the quantum Monte Carlo method suffers from the 
sign problem for frustrated systems. However, it is amazing that just owing to geomet- 
rical frustration some possibilities for rigorous analysis emerge. Recently, it has been 
recognized that many geometrically frustrated lattices (including the kagome lattice, the 
checkerboard lattice and the pyrochlore lattice) admit a simple class of exact eigenstates 
christened localized magnons [7|, |8j. These states become the ground states in strong 
magnetic fields and they are relevant for many physical properties of a wide class of 



lighly frustrated quantum antiferromagnets in the 



3, B S Q, H Q, Q. 
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ow-temperature strong-field regime 



22|, |23j . In particular, the localized- 



magnon states are responsible for magnetization jumps which the ground-state magnetiza- 



tion curve exhibits at the saturation field 
spin-Peierls lattice instabi 



saturation field 
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23], may lead to a high-field 



191 ]. and imply a residual ground-state entropy at the 
14J . Moreover, these states dominate the low-temperature thermo- 



dynamics in the vicinity of the saturation field [12 
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and may lead to an 



order-disorder phase transition of purely geometrical origin \U, 

This paper focusses on new developments concerning localized-magnon effects. We men- 
tion several reviews on this subject jj, 24, 25] which, however, do not cover some recent 



studies on the universal thermodynamic behavior which emerges at low temperatures around 
the saturation field. We will therefore concentrate on the results obtained during the last 



2 



two years 



21 



22). In passing, we will provide several complementary new results. 



Although the effects of the localized magnons have not been observed experimentally 



so far, recent studies on the spin-1/2 diamond-chain compound azurite Cu3(C0 3 ) 2 (0 



26. 



27| and the frustrated quasi- two-dimensional spin-1/2 antiferromagnet CS2CUCI4 



28] 



represent closely related physics. 

The present paper is organized as follows. In Section [Til we recall some basic results 
concerning the localized-magnon states which are then used throughout the paper. Further, 
in Section IIII} we discuss the linear independence of the localized-magnon states with the 
smallest localization area [21] and their completeness as a basis of the high-field ground 
states. Next, in Section IIVI we discuss the universal thermodynamic behavior which the 
considered frustrated quantum antiferromagnets exhibit at low temperatures in the vicin- 
ity of the saturation field 
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221 ] . We distinguish different types of the universal behavior 
depending on the specific classical lattice-gas model which represents the low-energy de- 
grees of freedom of the spin system. We focus on several representative models, namely, 
the diamond chain (hard-monomer universality class), the frustrated two- leg ladder and 
the kagome-like chain (one-dimensional hard-dimer universality class), and the frustrated 
bilayer lattice (two-dimensional hard-square universality class). Finally, in Section |V] we 
summarize and briefly comment on some possibilities of experimental observation of the 
localized-magnon effects. 



II. LOCALIZED-MAGNON STATES 



In this paper we study the Heisenberg model on several geometrically frustrated lattices. 
Some of these lattices are shown in Figs. [H [2j The Heisenberg Hamiltonian of iV quantum 
spins of length s reads 

H=J2 J nm (\ (S+S~ + S~S+) + A«J - h S* . (1) 
(nm) ^ ' 

Here the sum runs over the bonds (edges) which connect the neighboring sites (vertices) on 
the spin lattice (see Figs, [fl [2j), J nm > are the antiferromagnetic exchange integrals between 
the sites n and m, A > is the exchange interaction anisotropy parameter, h is the external 
magnetic field, and S z = J2 n s n * s ^ ne ^-component of the total spin (magnetization). Note 
that the operator S z commutes with the Hamiltonian ([T|). For the examples shown in 
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FIG. 1: Three examples of one-dimensional spin lattices admitting localized magnons: (a) the dia- 
mond chain (hard-monomer universality class), (b) the frustrated two- leg ladder (one-dimensional 



hard-dimer universality class), (c) the kagome-like chain, originally introduced in 29f] (one- 
dimensional hard-dimer universality class). 



Figs. [T^l, [T|d, and [2k the exchange integrals take two values, namely, J2 for the vertical 
bonds and J\ for all other bonds; for the other lattices shown in Figs. [U [2] they are uniform 
Jnm — J ■ The exchange integrals J\ and J2 may also be assumed to satisfy certain relations, 



see, e.g., Ref. 



20| and below. Although the present analysis can be performed for arbitrary 



values of s and A, in what follows we concentrate on the extreme quantum case s = 1/2 
and isotropic interactions A = 1, in particular while performing exact diagonalization for 
finite systems |30fl. 



From Refs. [3|, [7|, |8|, |25( we know that the Heisenberg antiferromagnet ([I]) on the lattices 
shown in Figs. (H [2], as well as on some other lattices like the dimer-plaquette chain, the 



sawtooth chain, another kagome-like chain 
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3, 



39( (one-dimensional systems) , the square- 



kagome lattice, the star lattice (two-dimensional systems), and the pyrochlore lattice (three- 



FIG. 2: Three examples of two-dimensional spin lattices admitting localized magnons: (a) the 
kagome lattice (hard- hexagon universality class), (b) the checkerboard lattice (large-hard-square 
universality class), (c) the frustrated bilayer lattice (hard-square universality class). We also show 
auxiliary lattice-gas models (hard hexagons on a triangular lattice, large hard squares on a square 
lattice and hard squares on a square lattice) which describe low-energy degrees of freedom of the 
spin models in strong magnetic fields. 
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dimensional system), support localized-magnon eigenstates. The magnon may be localized 
on the (vertical) bond as for the lattices shown in Figs. [T^i, [lb, and [2b, on the V-part of the 
sawtooth chain, or on the even polygon (hexagon, square etc.) from which other lattices 
are built (kagome, pyrochlore, checkerboard, square-kagome etc.), see Figs. [Tb, [2b., and [2b. 
Then the explicit expression for the localized-magnon state reads 

|llm) = |lm),|s,...,s) e , (2) 

where \lm)i = (l/v / 2)(|s)i|s — 1)2 — \s — 1)1)5)2) for the lattices in which the magnon is 
trapped on the vertical bond or |lm); oc Y^m=i(~^) m s m\ s ) 1 • • • \ S ) L ^ or ^ ne lattices in which 
the magnon is trapped on the L-site polygon and \s, . . . , s) e denotes the ferromagnetically 
polarized environment. By direct calculation one can convince oneself that the state (TSJ) 
is indeed an eigenstate of the Hamiltonian (TjQ) with the eigenvalue E'fm — £1 where Efm 
is the energy of ferromagnetically polarized lattice and e diamond = S (J 2 + A(2J X + J 2 )), 

e ladder = ^ + A(4Jj + ^agomc = ^ + ^ flayer = ^ + + ^ for 

the diamond chain, the frustrated two-leg ladder, the kagome-like chain, and the frustrated 
bilayer lattice, respectively (here we put h = in (pQ)). 

Alternatively, one may diagonalize the Hamiltonian (TjQ) in the one-magnon subspace with 
S z = Ns — 1 to find that one of the magnon excitation branches is flat (dispersionless). For 
the diamond chain, the frustrated two-leg ladder, the kagome-like chain, and the frustrated 
bilayer lattice the dispersionless one-magnon energy is given by A^ iamond = —s(J 2 + A (2 J\ + 
J 2 )) + h, A 1 ^ = -s{J 2 + A(4J a + J 2 )) + h, A^ gom6 = -2s(l + 2A) J + h, and A^ ilayer = 
— s(J 2 + A(8Ji + J 2)) + h, respectively (note the correspondence of A k with ei). For some 
lattices (the sawtooth chain, the kagome-like chains, the kagome lattice, the checkerboard 
lattice etc.) the dispersionless magnon band is the lowest one, for other lattices it may 
become the lowest one if certain relations on the antiferromagnetic exchange constants are 
imposed. In particular, for the diamond chain and the frustrated two-leg ladder we have to 
assume J 2 > 2J 1; whereas for the frustrated bilayer lattice we have to assume J 2 > 4J X . If 
equality in the imposed relations holds, the dispersive higher-energy magnon band touches 
the dispersionless lowest-energy band at some values of the wave-vector (this also occurs for 
the kagome-like chains, the kagome lattice, and the checkerboard lattice, but not for the 
sawtooth chain). 

We pass to the subspaces with total S z = Ns — 2, . . . , Ns — n max . Here n max is the number 
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of the isolated localized magnons, each occupying the smallest possible area, for the closest 
packing. n max depends on the lattice and equals N/3, N/4, N/6, and N/4 for the diamond 
chain, the frustrated two-leg ladder, the kagome-like chain, and the frustrated bilayer lattice, 
respectively. Evidently, we can construct eigenstates of the Hamiltonian ([T]) of the form 

|21m) = \lm) h \lm) h \s,...,s) e , |n max lm) = \lm) h . . . |lm)/ nm Js, . . . , s) e , (3) 

which in the zero-field case h = have the energies Efm — 2ei, . . . , E F m — n max ei, respec- 
tively, provided that the trapping cells of the localized magnons cannot be directly connected. 
Obviously, there are many other eigenstates in each of these subspaces. The special impor- 
tance of the localized-magnon states is due to the fact that they may become ground states 
(or at least low- lying states) in their respective subspaces. In Refs. [7, 9] it was proven under 
some quite general assumptions that the localized-magnon states are indeed lowest-energy 
states in the corresponding sectors of S z = Ns, . . . , Ns — n max . More precisely, if we denote 
the minimal energy within the subspace with S z = Ns — n by E m i n (n), the following inequal- 
ity holds E min (n) > (1 - n)E FM + nE mhl (l) = E FM - n(E FM - E min (l)) = E FM - ne x for all 
n — 0,1, ... , 2Ns and spin-s Heisenberg systems with sufficiently general coupling schemes 
. The energy of n localized-magnon states is given by the expression on the r.h.s. of 
this inequality (for n = 1, . . . , n max ) and hence we conclude that the localized magnons are 
lowest-energy states in the subspaces with S z = Ns — 1, ... , Ns — n max . 

We note that the localized-magnon states (T5]), ([3]) are highly degenerate states. Obviously, 
the localized magnons of smallest area can be placed on a lattice in many ways. Moreover, 
for some lattices there are other states which have the same energy as the localized magnons 
of smallest area. In the following discussion we will concentrate on the checkerboard lattice 



and refer to Refs. [13] , \24\ for related remarks on the kagome lattice. First of all we notice 
that for many lattices, we can construct localized magnons occupying a larger area in addi- 
tion to the localized magnon of smallest the example in Fig. [3^. It is important to 
note that such eigenstates can be viewed as linear combinations of the simpler eigenstates 
corresponding to the localized magnons of smallest area. Moreover, in two and higher di- 
mensions there may be additional topological effects. In particular, for some lattices one can 
construct nested objects (defect states), compare the example in Fig. [3b. Furthermore, for 



such lattices there may be additiona 



boundary conditions are imposed 



localized magnons with nontrivial winding if periodic 
24\ . We shall come back to this issue when discussing the 



FIG. 3: The checkerboard lattice, (a) A localized magnon which occupies a larger area than the 
smallest possible one. (b) A three-magnon state which is not a large-hard-square state because of 
the nested 'defect' state at the lower left side of the lattice; one magnon is localized on each of the 
three loops. 



degeneracy of the ground states in the subspaces S z = Ns—2, . . . , Ns—n ma ^ for finite systems 
using exact diagonalization. In particular, we will see that the localized magnons of smallest 
area linearly span the total ground state space in the sectors S z = Ns — 1, ... , Ns — n max 
for some lattices (e.g., the diamond chain and the frustrated two- leg ladder with J 2 > 2J l5 
the sawtooth chain, or the frustrated bilayer lattice with J 2 > 4Ji). 

One is thus led to the task of enumerating all many-localized-magnon states. Let us 
denote the number of possibilities to put n isolated localized magnons, each occupying the 
smallest possible area, on a lattice with N sites by gjv(n) > 1. Furthermore note that the 
energy of a state with n localized magnons reads 

E n (h) = E FM - hNs -n(ei - h) (4) 

in the presence of an external magnetic field. Thus, at the saturation field h = hi = ex 
we find that the energy E n (h) does not depend on the number of localized magnons n. 
Hence, the total ground-state degeneracy at h — hi is at least W = y^"^ n x Qn M provided 



the considered localized-magnon states are linearly independent (see Ref. 21( and the next 



section). As a result of this degeneracy one obtains the following universal properties: 1) the 
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ground-state magnetization curve exhibits the jump at hi between the values Ns — n mSLX and 
Ns; this jump is accompanied by a plateau at the foot of the jump; 2) since W scales 
exponentially with iV the ground-state entropy per site at hi does not vanish but remains 
finite, namely at least klnW/N > 0. Moreover, 3) the lattice is unstable with respect to a 
deformation, which preserves the symmetry required for the existence of localized magnons 
and lowers the magnetic energy linearly with respect to the displacement. As a result, a 
field-tuned structural instability may take place. We do not discuss these issues further; 
more details can be found, e.g., in Refs. 0, Instead we turn to the localized-magnon 
effects at finite but low temperatures and strong magnetic fields around the saturation field. 



III. LINEAR INDEPENDENCE AND COMPLETENESS 



As discussed in the last section, the localized-magnon states constitute a highly degenerate 
ground-state manifold at the saturation field h = h\. One may expect that these states lead 
to a dominant contribution to the low-temperature thermodynamics for magnetic fields in 
the vicinity of the saturation field. The contribution of the localized-magnon states to the 
partition function can be written in the form 

Z Xm {T, h, N) = Y^g N (n) exp f-^A J = exp (- EpM hNS J ]T g N {n) exp (^™) (5) 

n=0 V / V / n=0 

where Eq. (@J was used and fi = e\ — h = hi — h. However, several questions, briefly touched 
already in the previous section, have to be discussed before determining the degeneracy of 
the localized-magnon states gN(n) and analyzing the thermodynamic properties on the basis 
of Eq. §}. 

First, we have to clarify whether the set of localized magnons of smallest area is linearly 
independent in each sector of S z = Ns, . . . , Ns — n max . Only if this is the case, all localized 
magnons of smallest area contribute to the partition function of the spin system and we do 
not have to take care about the states which are their linear combinations (like the state 
shown in Fig.JHk). Bearing in mind the thermodynamic limit N — > oo we need at least linear 
independence in all sectors of S z except probably in a few (finite number) sectors of S z . 

Although we do not have a general theory of linear independence or dependence of the 
localized-magnon states of smallest area, a detailed and systematic analysis for a wide class 
of lattices hosting localized magnons was presented in Ref. 2l| including rigorous proofs for 



lattices in one and two dimensions. We discuss briefly these findings below following Ref. 21]. 
In this discussion we will refer to a localized magnon of smallest area simply by "localized 
magnon". It appears convenient to group the frustrated lattices into several classes and to 
examine linear independence for each class separately. Thus, we distinguish the following 
classes: the orthogonal type (the diamond chain (Fig. Wfa), the dimer-plaquette chain, the 
frustrated two- leg ladder (Fig. [lb), the square-kagome lattice, the frustrated bilayer lattice 
(Fig. [2b)); the isolated type (the sawtooth chain, the kagome-like chain I (Fig. [Tb) and the 
kagome-like chain II); the codimension 1 type (the kagome lattice (Fig. [2b,), the star lattice, 
the checkerboard lattice (Fig. Eh)); and the higher codimension type (the pyrochlore lattice). 

To examine the linear independence of a finite sequence of k vectors in some Hilbert 
space it is useful to introduce the k x k Gram matrix G = | \Gy\ | of all scalar products. The 
rank of G equals the dimension of the linear span of the considered set, i.e., the number of 
linearly independent states, and the set is linearly independent iff detG > 0. We call the 
dimension of the null space of G the codimension and hence the codimension will equal the 
number of independent linear relations between n localized- magnon states. 

It can be proven (theorem 1 of Ref. |2lJ]) that if the set of localized- magnon states is 
linearly independent in the sector S z = Ns — 1, then it is also linearly independent in 
the sectors S z = Ns — n, n = 2,...,n max . Therefore, in many cases it is sufficient to 
consider only localized one-magnon states. In particular, for the lattices which belong to 
the orthogonal type, the cells in which the localized magnons are concentrated are disjoint 
and any two different localized one-magnon states are orthogonal. As a result, Gij oc 5ij 
and according to the Gram criterion all localized one-magnon states (and hence all localized 
n-magnon states for every n = l,...,n max ) are linearly independent. Consider next the 
lattices which belong to the isolated type. It can be proven (theorem 3 of Ref. [21|) that if 
the trapping cell contains a site which is not contained in any other trapping cell ( "isolated 
site" ) then all localized one-magnon states (and hence all localized n-magnon states for every 
n — 1, . . . , n max ) are linearly independent. 

We pass to the lattices of the codimension 1 type. Suppose that every spin site is contained 
in exactly two different trapping cells (this is just the case of the lattices of the codimension 
1 type, e.g., the kagome lattice (Fig. [2b), the checkerboard lattice (Fig. [2b) etc.). Then it 
can be proven that there is at most one linear relation between localized one-magnon states 
(theorem 4 of Ref. 2l|). Moreover, in such a case the set of localized-magnon states in the 
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subspaces S z = Ns—n, n = 2, . . . , n max is linearly independent (theorem 5 of Ref. [211]) . Note 



that boundary conditions are crucial in this case. As was pointed out already in Ref. [24|, 
there is indeed one relation between localized one-magnon states for the kagome lattice 
if periodic boundary conditions are imposed. For the checkerboard lattice with periodic 
boundary conditions there is a similar linear relation 21|. By contrast, all localized one- 
magnon states are linearly independent if one considers the kagome or checkerboard lattice 
with open boundary conditions. Note that for the pyrochlore lattice there is more than one 
linear relation between one-magnon states localized on hexagons and the linear relations 
between such localized- magnon states exist also in the subspaces S z = Ns — n for n > 1 



2i 



The second important issue concerns the completeness of the localized-magnon states in 
each sector of S z = Ns, . . . , Ns — n max . More precisely, a quantitative description of the 
low-temperature thermodynamics of the spin model in the vicinity of the saturation field 
is obtained from the smallest-size localized magnons via Eq. ([5l) only if there are no fur- 
ther thermodynamically relevant contributions. Again we do not know a general answer. 
However, there are topological arguments supporting the existence of such additional con- 
tributions in certain higher- dimensional systems, and suggesting their absence otherwise. 
Furthermore, we can perform a quantitative comparison for each specific spin system sepa- 
rately with numerical results for finite systems. 

From exact diagonalization data for finite systems with periodic boundary conditions we 
infer that for some of the lattices there are no additional states, or at most one or two 
additional states in the subspace S z = Ns — 1. A perfect correspondence is found for the 
diamond chain and the frustrated two-leg ladder with J 2 > 2 J\ as well as the sawtooth chain 
and certain two-dimensional models with finite localization regions such as the frustrated 
bilayer lattice with J 2 > 4J X (see Ref. [2oJ and Tabled]). In other cases such as the kagome- 
like chains there are only one or two extra states which can be traced to one-magnon states 
from the dispersive band which at one or two values of the wave-vector have the same 
energy as the states from the flat band. In these cases which include in particular the one- 
dimensional models, a quantitatively accurate description is expected from the smallest-size 
localized magnons. 

In contrast, the topology of certain two-dimensional lattices permits the construction of 
additional, nested localized-magnon states [13l. l24|. An examples of such a 'defect' state is 
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TABLE I: The degeneracies and the energy gaps for various finite s = 1/2 frustrated bilayer lattices: 
exact diagonalization data for finite systems with J2 = 5 J\ vs hard-square predictions. TV" is the 
number of sites in the spin lattice, DGS is the degeneracy of the ground state as it follows from 
exact diagonalization for a given S z , # HSS is the number of configurations with N/2 — S z hard 
squares, A is the energy gap between the ground state and the first excited one in units of J\. The 
dots for N = 64 indicate omitted sectors with 16 < S z < 26. 
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shown in Fig. [3b. Such defect states occur first in the two-magnon sector and, once again, 
the issue of linear independence from the smallest-area localized magnon states arises. We 
would like to point out that boundary conditions again play a crucial role: for open boundary 
conditions, the two-magnon state obtained from the two nested loops at the lower left corner 
of Fig. [3b is indeed a new state. However, for the periodic boundary conditions which we 
are considering here, such states can be expressed as a linear combination of smallest-area 
localized two-magnon states. Let us explain this in the case of the checkerboard lattice. 
First note that in a finite checkerboard lattice with periodic boundary conditions, a magnon 
state localized on a closed loop L such as the large loop in Fig. [3b can be written as a linear 
combination of all magnon states localized on the smallest squares contained in L. But this 
linear combination is not unique since the set of all smallest localized one-magnon states 
is not linearly independent. Another linear combination yielding the same state could run 
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TABLE II: The degeneracies and the energy gaps for the s = 1/2 kagome lattice: exact diagonal- 
ization data for finite systems vs hard-hexagon predictions. N is the number of sites in the spin 
lattice, DGS is the degeneracy of the ground state as it follows from exact diagonalization for a 
given S z , # HHS is the number of configurations with N/2 — S z hard hexagons, A is the energy 
gap between the ground state and the first excited one in units of J. Dots indicate some sectors 
with S z > 7 N/18 which have been omitted for larger values of N. For the two-magnon sector one 
has # HHS = iV 2 /18 - 7 JV/6 and DGS = iV 2 /18 - N/2 + 1. 
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over all magnon states \q) localized on the smallest squares q not contained in L. Consider 
in particular the octagonal loop L containing a smallest square go with a magnon state |go) 
in its center shown in Fig. ^jp. Then the state of the two magnons localized on L and go can 
be written as a linear combination of large-hard-square two-magnon states |g) |go) where |g) 
runs over all one-magnon states of smallest area not contained in L. 

Note that the derivation of the linear relation relies on the absence of obstacles such 
as open boundaries or other localized magnons outside the two-magnon defect state. We 
therefore believe that the majority of many-magnon states constructed with these defect 
states (in particular the three-magnon state of Fig. ^jp) is linearly independent of many- 
magnon states containing only those of smallest area. The many-defect states are therefore 
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TABLE III: The degeneracies and the energy gaps for various finite s = 1/2 checkerboard lattices: 
exact diagonalization data for finite systems vs large-hard-square predictions. N is the number 
of sites in the spin lattice, DGS is the degeneracy of the ground state as it follows from exact 
diagonalization for a given S z , jf= LHSS is the number of configurations with N/2 — S z large-hard- 
squares, A is the energy gap between the ground state and the first excited one in units of J. Dots 
indicate some sectors with S z > 3N/8 which have been omitted for larger values of N. For the 
two-magnon sector one has # LHSS = N 2 /8 - 9 N/4 and DGS = N 2 /8 - 5 N/4 + 1. 
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expected to yield another finite (even if small) contribution to the ground-state entropy at 
the saturation field. Indeed, numerical data for finite lattices exhibits a larger ground-state 
degeneracy than predicted from the effective hard-object description (see Table [Til for the 
kagome lattice and Table [TTTI for the checkerboard lattice). This difference should remain 
relevant in the thermodynamic limit according to the preceding argument based on the 
defect states. 

Before we proceed, we would like to comment on some specific sectors for the kagome lat- 
tice (Table HTj) and the checkerboard lattice (Table [TTTI) . Firstly, the ground-state degeneracy 
(DGS) in the one-magnon sector is A^/3 + l (N/2 + 1) for the kagome (checkerboard) lattice. 
There are two possible interpretations (see also 24|): in momentum space this corresponds 
to the N/3 {N/2) states of the flat branch plus one additional state where the next disper- 
sive branch touches the flat branch; in real space this corresponds to iV/3 (N/2) smallest 
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localized-magnon states subject to one linear relation plus two additional states which wind 
once around the boundaries. Secondly, the exact degeneracy in the two-magnon sector is 
given by iV 2 /18 — N/2 + 1 for the kagome lattice and N 2 /8 — 5 iV/4 + 1 for the checkerboard 
lattice. The difference with respect to the number of smallest-area localized-magnon config- 
urations is again due to states with non-trivial winding. Let us explain this briefly for the 
checkerboard lattice. First, there are N configurations of the form \q)\w x ) or \q)\w y ) where 
q runs over all squares and w a is an arbitrary loop which winds around boundary a = x, 
y sufficiently far away from q. As we will explain in detail elsewhere, taking into account 
states \w a )\w' a ) with double winding together with one special diagonal state and the linear 
relations between these states, one finds N + 1 additional linearly independent two-magnon 
states with non-trivial winding, i.e., exactly the same number as observed numerically. The 
difference 2 A^/3 + l between the number of configurations of two hard hexagons and the exact 
ground-state degeneracy in the two-magnon sector for the kagome lattice can be explained 
in an analogous way. Finally, let us look at the sector with the closest packing of localized 
magnons, i.e., S z = 7 N/18 for the kagome lattice and S z = 3 N/8 for the checkerboard 
lattice. For the N = 40 checkerboard lattice we find the same number of ground states and 
large-hard-square states for S z = 15, as expected for a closest packing. By sharp contrast, 
the kagome lattice gives rise to 8 ground states for N = 36, and 4 for N = 45, 54 in the 



sector with S z = 7 N/18 while one expects only 3 for the magnon crystal [3J, |§J, |24|. Since 
the N = 36, 45 and 54 kagome lattices should be sufficiently large to eliminate boundary 
artifacts, the origin of the additional state(s) is unclear at present. In particular, it remains 
to be clarified whether (essentially) all ground states are described by localized magnons if 
all topological non-trivial configurations (including defect states and states with non-trivial 
winding) are properly accounted for. 

One further issue is whether the ground-state manifold is separated from other states by 
a finite energy gap. We can again draw some conclusions concerning the energy gap from 
exact diagonalization data. In previous papers 20, |22j we introduced a measure for the 
thermodynamically relevant energy separation Ados between the ground-state manifold and 
the other eigenstates of the system. For simplicity, in the present paper we report the energy 
gap A between the ground-state energy and the next smallest energy level in each sector 
S z . Tables H] HT] and IIHI present the values of A for some finite frustrated bilayer, kagome, 
and checkerboard lattices, respectively. For the frustrated bilayer lattice with J 2 = 5 J\, 
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the localized magnons are separated from all higher excitations by an energy gap A > J\ 
(see Table HJ. For the kagome lattice one has the additional complication that there is no 
gap to the next branch of excitations in the one-magnon sector. Nevertheless, two-magnon 



scattering states were estimated to have an energy gap A w 0.24 J [13|, [24|. While this may 
be a valid estimate in the two-magnon sector of the kagome lattice, in higher sectors there 
are definitely excitations at substantially lower energies (see Table HTjl . In fact, analysis of 
further excited states (not shown here) indicates the onset of a thermodynamically relevant 
density of states at energies of the order of only 1CT 2 J. 



IV. LOW-TEMPERATURE STRONG-FIELD THERMODYNAMICS. LATTICE- 
GAS DESCRIPTION 

As mentioned above, the smallest-area localized-magnon states may dominate the low- 
temperature thermodynamics in the vicinity of the saturation field. After having checked 
their linear independence we would like to discuss their contribution to the canonical parti- 
tion function of the spin system in more detail. We start from Eq. (jSj), where this contribution 
is given. We emphasize once again that this formula describes the low-temperature thermo- 
dynamics near the saturation field accurately provided that (i) there are no other ground 
states (apart from the smallest-area localized-magnon states) in the corresponding sectors 
of S z or that the contribution of such extra states is vanishingly small as N — > oo, and 
that (ii) excited states in these sectors are separated by a finite energy gap from the ground 
states. In Eq. (jSJ) gN( n ) is the degeneracy of n isolated smallest-area localized magnons on 
a spin lattice of N sites. It is useful to consider #Ar(n) as the canonical partition function 
Z(n, AT) of n hard-core objects on an auxiliary lattice of Af oc N sites (Af = N/3, N/2, N/3, 
and N/2 for the diamond chain, the frustrated two-leg ladder, the kagome-like chain, and 
the frustrated bilayer lattice, respectively). We can write the grand-canonical partition func- 
tion of hard-core objects on this lattice as E(T,[i,Af) = J^=q x 9n(ji) exp (fin/kT), where 
/i is the chemical potential of the hard-core objects. The simple reason why a hard-core 
object lattice-gas description emerges here is the existence of the "hard-core rules" which 
the localized-magnon states must respect in order to be eigenstates of the spin Hamilto- 
nian. Note that these rules differ for various spin lattices: the diamond chain of Fig. [TK 
gives rise to hard monomers, the frustrated ladder (Fig. [lb) and the kagome-like chain 
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(Fig. [lb) are described by hard dimers, while the hard-core objects for two-dimensional lat- 
tices are illustrated in Fig. [2j Using S(T, fi,Af) we arrive at the basic relation between the 
localized-magnon contribution to the canonical partition function of the spin model and the 
grand-canonical partition function of the corresponding hard-core object lattice-gas model, 

Z\ m (T, h, N) = exp (- EFM ~ T hN ^ E{T,n,AT) , (6) 

with jj, — h\ — h. From Eq. ([HD we find the Helmholtz free energy of the spin system 
F lm (T,h,N) E FM M In S(T,/i, A/") 

77- = — 77- hs — kl — —. . 7 

N N N M w 

The entropy S, the specific heat C, the magnetization M = (S z ), and the susceptibil- 
ity x follow from (J7J) according to usual formulae, S\ m (T, h, N) = —dF\ m (T,h,N)/dT, 
C lm (T,h,N) = TdS lm (T,h,N)/dT } M lm (T, h, N) = Ns - kTd\nE(T,fi,Af)/d^, and 
Xhn{T, h,N) = dMi m (T, h, N)/dh, respectively. Note that the thermodynamic quantities 
of the spin system depend on the temperature T and the magnetic field h essentially only 
through the combination x = (hi — h)/kT j^. 

To test the hard-object description, we have performed full diagonalization of spin-1/2 
isotropic Heisenberg systems (i.e., A = 1 in Eq. (DO)), imposing periodic boundary conditions. 
For the diamond chain and the frustrated ladder we have also exploited the local conservation 
laws. First we rewrite the Hamiltonian in terms of the total spin on the vertical dimers in 
Fig. Q] 32|. For a local spin 1/2, each dimer can only be in the singlet or triplet state. 
A singlet cuts the system into smaller fragments. It is therefore sufficient to compute the 
spectra of one periodic fragment where all vertical dimers are in the triplet state, and smaller 
open fragments where all consecutive dimers are again in the triplet state. In this manner 
it requires only a moderate effort to obtain the complete spectra for a diamond chain and a 
frustrated ladder with = 24, while A^ = 24 would be inaccessible for the diamond chain 
with a full diagonalization of the original model. 



A. Hard-monomer universality class 

First we consider the hard-monomer universality class which includes the diamond chain, 
the dimer-plaquette chain and the square-kagome lattice. The hard-monomer restriction for 
these lattices means that it is forbidden to have two (or more) localized magnons in the 
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same trap. We focus on the diamond chain (see Fig. Wfa), and refer the interested reader 



to Ref. 



20| | for the other lattices. A straightforward calculation yields the grand-canonical 



partition function of a gas of hard monomers 

~CZ>,A0= (l+exp^f. (8) 

Explicit analytic expressions for thermodynamic quantities can be obtained easily [^(J . Note 
that the thermodynamic quantities per site are independent of the system size N. At the sat- 
uration field we have a residual ground-state entropy S\ m (T, hi, N)/kN = (Af/N) In 2. The 
specific heat C\ m (T, h l} N)/kN exhibits two identical maxima of height rs 0.43922884A/" /N 
at x « ±2.39935728. 

Some typical dependencies of the thermodynamic quantities on the field and the temper- 
ature for the diamond chain with Ji = 1, J 2 = 3 and the corresponding hard- monomer data 
obtained on the basis of Eq. (jBJ) are shown in Fig. HI left column. Deviations between exact 
diagonalization (ED) and hard monomers (HM) are observed only in the specific heat C, 
and only for kT > 0.2 J\ (lowest panel in the left column of Fig. H]). Note furthermore that 
thermodynamic quantities are symmetric under x — > —x within the hard-monomer picture. 
In particular, the hard-monomer description yields a susceptibility x an d a specific heat C 
for the diamond chain which coincide at h = 3.8 J\ and h = 4.2 J]_. 



B. One-dimensional hard-dimer universality class 

The frustrated two-leg ladder, the kagome-like chains, and the sawtooth chain belong to 
the one- dimensional hard-dimer universality class, i.e., the rules for the localized magnons 
obey the restrictions for rigid dimers on a one- dimensional lattice: each trapping cell can 
only be occupied by one localized magnon and neighboring trapping cells cannot be simul- 
taneously occupied by localized magnons. We focus on the frustrated two-leg ladder and 
the kagome-like chain of type I (see Figs, [lb, Hfc). The grand-canonical partition function of 
one-dimensional hard dimers can be calculated with the help of the transfer-matrix method: 

2(T,/,,A/-) = Af + Af, *i.2 = i±^ + expA (9) 

Explicit analytic expressions for thermodynamic quantities can again be obtained easily. 
Note that for hard dimers the thermodynamic quantities per site depend on the size M. In 
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FIG. 4: M(T, h, N) /N, X (T, h, N) /N, and S(T, h, N) /kN vs h at low temperatures; x(T, h, N) /N 
and C(T,h,N)/kN vs kT around the saturation field. From left to right: diamond chain with 
J\ = 1, J2 = 3, frustrated two-leg ladder with J\ = 1, J2 = 3, kagome-like chain with J = 1. 
We set the field range foi — 1 < /i < /ii + 1 and the temperature range < kT < 0.5. The exact 
diagonalization (ED) data (symbols) refer to finite systems of sizes N = 24 (diamond chain and 
frustrated ladder) and N = 18 (kagome-like chain). The analytical predictions for hard monomers 
(HM) and one-dimensional hard dimers (HD) (N — > 00) are shown by lines. 
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the thermodynamic limit only the 



this limit can be found in Refs. 
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are 



icit expressions in 



20] for the general 



rgest eigenvalue Ai plays a role. Exp 
24| for the sawtooth chain and in Ref. 
case. At the saturation field one finds a residual ground-state entropy Si m (T, hi, N)/kN = 
{N/N) ln((l + v / 5)/2). The specific heat C lm (T, hi, N)/kN exhibits two maxima of height 
« 0.34394234 N/N (at x « -2.81588498) and » 0.26887020 Af/iV (at x « 4.05258891). 
Some typical dependencies of the thermodynamic quantities on field and temperature are 
shown in Fig. HI middle and right columns for the frustrated ladder with J\ = 1, J 2 = 3 and 
the kagome-like chain with J — 1, respectively. Note that the hard-dimer results are not 
symmetric around the saturation field, i.e., all thermodynamic quantities are different for x 
and —x. 

For the frustrated ladder with J?, = 3 J\, we observe again systematic differences between 
exact diagonalization and hard dimers in the specific heat C at high temperatures (lowest 
panel in the middle column of Fig. H]). The remaining differences for low temperatures and 
h < h\ are due to finite-size effects, since in Fig. H] we compare finite spin systems with 
infinite hard-dimer systems. If the comparison is performed for the same system size j^oj], 
better agreement can be observed. One important source of finite-size effects is a two-fold 
degeneracy of the ground state at M/N =1/4. This ground-state degeneracy is evident in 
the finite value of the ED results for the entropy S at the left side of the second panel in 
the middle column of Fig. HI 

In the case of the kagome-like chain (right column of Fig. HJ) we find good agreement 
between exact diagonalization and hard dimers on the high-field side h > hi and sufficiently 
low temperatures. However, on the low-field side h<hiwe observe even stronger deviations 
than for the frustrated ladder. Indeed, there are two steps at h 2 ~ 2.67 J and h 3 « 2.11 J 
in the zero-temperature magnetization curve of the iV = 18 kagome chain with M/N < 1/3 
(see also Fig. 2 of Ref. {23]). Since the region M/N < 1/3 cannot be described with hard 
dimers, the region with h < 2.7 J falls completely outside the validity of the hard-dimer 
picture. 

C. Two-dimensional lattice gases, hard-square universality class 

The frustrated bilayer lattice, the kagome lattice as well as the checkerboard lattice are 
even more interesting, since the corresponding two-dimensional classical hard-core systems 
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exhibit a second-order finite-temperature order-disorder phase transition [40l. |41|. |42||. How- 
ever, recall from Section HTT1 that the spin model exhibits extra ground states on the kagome 
and checkerboard lattices which cannot be described by hard hexagons or large-hard-squares, 
respectively. Moreover, the gap to excited states is pretty small, see Tables HT1 and IIHI There- 
fore, the description of the kagome and checkerboard lattices in terms of hard-core objects 
is expected to be only a qualitative one. 

We therefore focus on the frustrated bilayer antiferromagnet as an example for a finite- 
temperature order-disorder phase transition 22j . In this case, exact diagonalization provides 
clear evidence that all ground states are mapped onto the hard-square configurations on a 
square lattice, see Table fl] Moreover, the gap to excited states is large. Thus, the low- 
temperature strong-field thermodynamics should be determined completely by the hard- 
square problem. Although we do not know the exact analytical result for the grand thermo- 
dynamical potential — kT In SfT, wJV") /M of hard squares on a square lattice, the properties 



of the model are well known 
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42] . In particular, the hard-square model exhibits a phase 



transition at z c = exp(/i c /A;T) = 3.7962 . . . between the low-density phase [z < z c ), in which 
both sublattices of the underlying square lattice are equally occupied, and the high-density 
phase (z > z c ), in which one of the sublattices becomes more occupied than the other one. In 
spin language, the phase transition has a purely geometrical origin and indicates the order- 
ing of localized-magnon states as their density varies with field or temperature. The phase 
transition belongs to the two-dimensional Ising universality class. Hence, the specific heat 
should show a logarithmic singularity at the critical point. Fig. [5] shows results for the tem- 
perature dependence of the specific heat around the saturation field for the frustrated bilayer 
lattice. The available exact diagonalization data are restricted to rather small spin systems. 
However, they demonstrate a perfect agreement with the results for the corresponding finite 
hard-square model: for J 2 = 5 Ji and h = 0.99 hi = 8.91 Ji or h = 1.01 h% = 9.09 Ji, both 
data sets in Fig. [5] are indistinguishable for temperatures kT < 0- 1 J\- Bigger hard-square 
systems can be studied using classical Monte Carlo simulations 43|. The left panel of Fig. [5] 
shows that a logarithmic singularity in the dependence of C vs T develops on the low-field 
side with increasing size of the hard-square model. Since this transition occurs within the 
temperature region where the exact diagonalization data are perfectly reproduced by the 
hard-square model, we expect that this singularity also appears in the spin model in the 
thermodynamic limit iV — > oo. 
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FIG. 5: C(T,h, N)/kN vs kT around the saturation field h\ = 9 for the frustrated bilayer lattice 
with J\ = 1, J2 = 5. The exact diagonalization data (filled symbols) refer to a finite spin system 
of ./V = 20 sites. The analytical results (empty symbols) refer to a finite hard-square system of 
N = 10 sites. The Monte Carlo simulation data (lines) are obtained for the hard-square system 
on finite lattices with M up to 800 x 800. 



To summarize this subsection, the frustrated bilayer lattice provides an example where 
a two-dimensional hard-core lattice gas completely covers all low-energy states of the spin 
model. In this model, there is a clear phase transition which corresponds to a crystallization 
of hard squares, i.e., magnons localized on the vertical dimers. We would like to emphasize 



44j since only 



that such a phase transition does not contradict the Mermin- Wagner theorem 
a discrete symmetry is broken spontaneously. This demonstration of a finite-temperature 
phase transition in a two-dimensional interacting many-body spin model is an interesting 
example for the impact of the localized-magnon states on the physical properties of a wide 
class of frustrated magnets. 



D. Region of validity 



Finally, we add some general remarks about the region of validity of the hard-core lattice- 
j;as description. Obviously, this effective picture of the spin model is accurate only in some 
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region in the h — T-plane around the point h = hi, T = 0. The field h < hi, until which 
the hard-core object picture should work at T = 0, is related to the width of the plateau 
hi — h 2 preceding the jump in the ground-state magnetization curve (si, 10, 11, 23] with h 2 
being the difference between the ground-state energy in the sectors S z = Ns — n max and 
S z = Ns — n max — 1. We could also try to estimate a characteristic temperature T*{h), below 
which the hard-core object picture should work at a certain magnetic field h. At h = hi we 
can find T* from the temperature dependence of the specific heat C. Indeed, the hard-core 
object prediction is C = for all temperatures at h = hi. However, as can be seen in the 
corresponding panels in Fig. H]this is a valid approximation for kT < 0.2Ji (kT < 0.15Ji) 
for the diamond chain (frustrated ladder), whereas the corresponding temperature region 
is much smaller for the kagome-like chain. For h > hi the temperature T* depends on the 
difference h — h\. Numerical data for finite systems indicates that T* increases with growing 
h — h\. 

Lastly, we note that a description of the low-energy degrees of the spin systems can 
be extended by relaxing the hard-core rules, e.g. by rendering the infinite nearest-neighbor 



repulsion finite or permitting double occupation of the auxiliary lattice sites [20|, |24j . In such 
a case we may achieve a better agreement with exact diagonalization data in a wider range 
of parameters around the point h = hi, T = 0, loosing, however, the universal dependence 
on h and T only via the parameter x = (hi — h)/kT. 



V. CONCLUSIONS 



In this paper, we have discussed the universal properties of some highly frustrated quan- 
tum Heisenberg antiferromagnets supporting localized-magnon eigenstates. Universal be- 
havior emerges owing to the localized-magnon states which become the ground states around 
the saturation field and can be separated from the higher-energy states by an energy gap. 
We find several universality classes depending on the specific lattice-gas model of hard-core 
objects which describes the low-energy degrees of freedom of the spin model in strong mag- 
netic fields. For the one-dimensional models, the lattice gas yields a quantitative description 
of the thermodynamics of the full spin model close to the saturation field and at sufficiently 
low temperatures. 

Higher dimensions may be even more interesting since they allow for a finite-temperature 
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crystallization phase transition of the hard-core objects. As a two-dimensional example, we 
have focused on the hard-square universality class which contains, e.g., the frustrated bilayer 
quantum Heisenberg antiferromagnet in the region with sufficiently strong interlayer coupling 



221 ] . Such a spin model exhibits an order-disorder phase transition of a purely geometrical 



origin which reflects the geometrical ordering of the localized magnons. It turns out that the 
phase transition of hard squares belongs to the two-dimensional Ising universality class and 
is characterized by a logarithmic singularity of the specific heat just below the saturation 
field. 

New numerical results and complementary arguments presented in this paper indicate 
that the situation is more complicated in the case of the kagome and checkerboard lattices. 
Here, there seem to be thermodynamically relevant contributions to the ground-state man- 
ifold beyond that of hard hexagons and large-hard squares, respectively. Accordingly, in 
these cases the gas of smallest hard objects cannot be expected to yield a quantitatively ac- 
curate description of the spin model in any parameter regime. Still, the universality class of 
a possible crystallization phase transition should be predicted correctly by such an effective 
low-energy theory. 

We believe that the properties of localized magnons elucidate the physics of frustrated 
quantum antiferromagnets in high magnetic fields and thus are useful for a general under- 
standing of related compounds. Even more, recent experiments on the spin-1/2 (distorted) 
diamond-chain compound azurite Cu^fCO^WOH )^ j^M Q and the frustrated quasi-two- 
dimensional spin-1/2 antiferromagnet CS2CUCI4 [28| raise hopes for a direct comparison 
with the theoretical models discussed in this paper, although so far there is no clear ex- 
perimental observation of the pronounced quantum effects predicted for low-dimensional 
spin-1/2 antiferromagnets with localized magnons yet. 

The Cu 2+ ions in azurite Cu^CO.sWOH^ form infinite chains with the structure of a 
spin-1/2 distorted diamond chain 26|, |27j. The high-field magnetization M(h) of azurite 
exhibits a plateau at 1/3 of the saturation magnetization and a further steep increase as 
the magnetic field tends to the saturation value of about 32.5 T 26| , |27| . Note that an 
experimentally accessible saturation field is an attractive feature of this compound. Fits 
of the magnetization curve and thermodynamic properties with high-temperature series 



371 ] and numerical results yield the following estimates for the exchange interactions |27| : 



J™ = 8.6 K, J( = 19 K, J2 = 24 K (in Fig.[TK the bonds with weaker interaction J™ run from 
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south-west to north-east whereas the bonds with stronger interaction J( run from north- 
west to south-east). Deviations from the ideal diamond chain geometry are not necessarily 
a major problem (see below). However, the quoted estimates for the exchange interactions 
are not in the region J 2 > 2 Ji which is required to render the localized magnons low- 
energy excitations (compare also the phase diagram of the distorted diamond chain 45]). 
Nevertheless, the values of the exchange couplings are still under debate [27]. 

CS2CUCI4 is a frustrated quasi-two-dimensional spin- 1/2 antiferromagnet with a low 
saturation field of only about 8.5 T Q]. For this compound, measurements of the low- 



d have already been 
. The dominant ex- 



temperature behavior of the specific heat around the saturation fie. 
performed, exhibiting a strong dependence on the magnetic field 
change interactions in CS2CUCI4 correspond to an anisotropic triangular lattice with strong 
interactions along one 'chain' direction. For such a model one can construct magnons local- 
ized on the strongly coupled chains, much in the same way as for the frustrated square lattice 
3|, l8|]. Due to the quantization of momenta transverse to the chain direction, these localized 
magnon states turn out to be high-field ground states for up to at least six coupled chains in 
the parameter regime relevant to CS2CUCI4 46| . The localized magnons cease to be ground 
states as one approaches the thermodynamic limit, but they remain low-energy excitations. 
However, even in a case such as the frustrated square lattice where magnons localized on lines 
are high- field ground states for all finite systems, these are not thermodynamically relevant, 
but instead magnetic order occurs below the saturation field 47J. Indeed, CS2CUCI4 exhibits 
a finite-temperature magnetic ordering transition below the saturation field 46| such that 
inter-plane coupling will have to be taken into account for a quantitative description of the 
low-temperature specific heat 28j . Nevertheless, the existence of localized magnons and the 
strong field-dependence of the specific heat of CS2CUCI4 [21J are both related to the strong 
frustration of the anisotropic triangular lattice. 

With respect to experiments such as those on azurite it is desirable to be not restricted 
to the "ideal geometry" allowing existence of exact localized-magnon ground states and to 
examine the "stability" of our results against small deviations from the relation which we 
have imposed on the exchange interactions. Numerical studies for finite spin systems [14J] 
suggest that the main features coming from the localized-magnon states survive in case of 
small deviations from the ideal lattice geometry. In general, we can argue that due to small 
deviations from the ideal geometry, the flat magnon band becomes slightly dispersive but 



25 



i a quantum 



481 ] . see also 



the hard-core constraint is preserved. Apparently, in this case we are faced wit 
hard-core object model (e.g., the quantum hard-square model studied in Ref. 
references therein). A study of the corresponding low-energy theories is beyond the scope 
of the present paper. 
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Note added: Two closely related preprints were submitted shortly after the 

present one [5l|. Ref. 49J confirms some of the new conclusions of our Section [TTTl Fur- 



thermore, Refs. 



49 



150] both contain Monte-Carlo results for the phase transition in the 
large-hard-square lattice gas. They agree on the location of the critical point while the 
universality class of the phase transition remains under debate. 
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